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Abstract 


Autonomous  amphibious  vehicles  have  the  potential  to  revolutionize  supply  chain  operations 
involving  the  transport  of  cargo  from  a  Sea  Base  to  troops  ashore.  In  2007,  a  team  of  Center 
for  Innovation  in  Ship  Design  (CISD)  interns  developed  a  design  concept  called  DUKW-21, 
an  amphibious  vehicle  meant  to  support  one  crew  member  onboard,  with  the  potential  for 
future  integration  of  an  autonomous  control  system.  While  DUKW-21  is  designed  to  be  a 
tracked  vehicle,  for  simplicity,  a  wheeled  1  /7th  scale  model  was  designed  and  built,  called 
DUKW-ling,  which  has  successfully  demonstrated  limited  autonomy. 

Unlike  an  air-cushion  vehicle,  a  tracked  or  wheeled  amphibious  vehicle  must  travel  through 
three  distinct  environments:  sea,  (dry)  land,  and  the  transition  zone,  where  the  vehicle's 
ground  propulsors  are  engaged  with  the  sea-floor,  but  the  vehicle  is  still  partially  submerged 
underwater.  In  this  report,  an  original,  baseline  model  of  the  dynamics  of  a  tracked  or 
wheeled  amphibious  vehicle  in  the  transition  is  developed.  This  model  is  then  used  to  study 
vehicle  routing  in  the  transition  with  the  objective  of  minimizing  travel  time.  In  particular,  a 
closed-form  characterization  of  a  preferred  travel  route  is  derived.  Under  certain  conditions,  it 
is  also  shown  that  the  baseline  model  of  the  transition  can  be  applied  to  hills  on  land. 
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Introduction 

Background 

The  United  States  military  currently  relies  heavily  on  large  transport  vessels  part 
of  the  Maritime  Pre-positioning  Force  (MPF),  and  developed  deep  water  port 
facilities  to  deliver  logistic  material  abroad.  The  difficulties  in  finding  and  securing 
port  facilities  in  hostile  areas  spurred  the  development  of  the  Sea  Base  concept; 
the  provision  of  mobile  port  facilities  in  controlled  waters.  The  key  challenge  to 
the  concept  is  the  development  of  a  supply  chain  to  transport  cargo  from  the  Sea 
Base  to  shore.  One  approach  is  to  use  a  medium  size  container  ship  to  transport 
cargo  from  the  Sea  Base  to  a  location  near  the  shore,  at  which  point  amphibious 
vehicles  can  run  continuous  delivery  missions  to  a  target  location  inland. 

Automating  the  amphibious  vehicles  in  such  an  application  would  have  significant 
advantages.  Flaving  the  capability  for  a  team  of  unmanned  amphibious  vehicles 
to  deliver  cargo  would  free  up  the  service  crew  for  other  tasks.  In  addition, 
unmanned  vehicles  can  maneuver  in  more  agile  ways  unconstrained  by  the 
acceleration  limits  that  humans  can  tolerate,  which  would  greatly  increase  the 
efficiency  of  the  cargo  transport  supply  chain.  The  Center  for  Innovation  in  Ship 
Design  (CISD),  a  part  of  the  Naval  Surface  Warfare  Center  Carderock  Division, 
has  been  developing  a  concept  for  such  an  amphibious  vehicle,  called  DUKW- 
21,  shown  in  Figure  1. 

In  2007,  a  team  of  CISD  interns  developed  a  design  concept  for  DUKW-21 ,  a 
tracked  amphibious  vehicle  with  a  human  support  system  for  a  single  crew 
member  that  would  facilitate  the  eventual  implementation  of  an  autonomous 
control  system  (Gonzalez  et.  al.,  2007).  Research  into  automating  DUKW-21 
started  in  2009,  when  a  CISD  intern  investigated  existing  ground  and  sea 
navigation  algorithms  and  ways  of  integrating  the  different  operational  modes 
(Flom,  2009). 

While  DUKW-21  is  designed  to  be  a  tracked  vehicle,  for  simplicity,  a  wheeled 
1  /7th  scale  model  was  developed,  called  DUKW-ling,  which  has  successfully 
demonstrated  limited  autonomy  (FAU,  2010).  In  this  report,  several  issues 
associated  with  a  tracked  or  wheeled  amphibious  vehicle  driving  up  the  sea-floor 
during  the  beaching  phase  are  studied  with  the  idea  that  in  the  future,  the  tools 
developed  here  can  be  tested  on  DUKW-ling. 
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Figure  1:  CISD  DUKW-21  concept,  borrowed  from  (Gonzalez  et.  al.,  2007) 

Challenges 

Unlike  an  air-cushion  vehicle,  a  wheeled/tracked  amphibious  vehicle  must  travel 
through  three  distinct  environments:  sea,  (dry)  land,  and  the  transition  zone, 
where  in  the  transition,  the  vehicle's  ground  propulsors  are  engaged  with  the  sea¬ 
floor,  while  the  vehicle  is  still  partially  submerged  in  the  water,  as  shown  in  Figure 
2.  We  call  the  border  between  the  sea  and  transition  the  engage  line  and  the 
border  between  the  transition  and  land  the  shoreline.  All  these  terms  are  defined 
more  precisely  later. 

Aside  from  additional  hydrodynamic  forces  acting  on  the  vehicle,  there  are  a  few 
characteristics  that  make  the  transition  distinct  from  dry  land.  Firstly,  the  tire/track 
to  ground  surface  friction  coefficient  is  generally  much  smaller  on  sand 
submerged  underwater  than  on  dry  land,  which  leads  to  a  significantly  smaller 
propulsive  force  that  the  vehicle  can  generate  in  the  transition.  This  makes 
optimizing  the  route  the  vehicle  follows  critical  because  poor  route  planning  can 
result  in  the  vehicle  being  immobilized  in  wet  sand,  which  would  increase  the 
time  it  takes  for  the  vehicle  to  complete  its  mission,  thereby  delaying  delivery  of 
the  supplies  that  the  vehicle  is  transporting. 

Another  factor  which  makes  the  transition  distinct  from  land  is  that  the  effective 
weight  of  the  vehicle  is  not  constant  due  to  the  change  in  buoyancy  as  the 
vehicle  surfaces.  This  implies  that  existing  navigation  algorithms  for  ground 
vehicles  are  not  necessarily  applicable  because  they  usually  assume  that  the 
weight  of  the  vehicle  is  constant. 


Figure  2:  Amphibious  operating  environments 
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Report  Overview 

We  make  several  contributions  in  this  report.  In  the  “Physical  Model  of  Transition” 
chapter,  we  propose  a  baseline  model  for  the  transition.  Then  in  the 
“Implications”  chapter,  the  model  is  used  to  study  vehicle  routing,  also  known  as 
trajectory  or  path-planning.  It  is  also  shown  that  under  certain  conditions,  the 
proposed  model  can  be  used  for  navigating  general  hills  on  dry  land.  This 
extends  the  application  of  our  research  from  autonomous  amphibious  vehicles  to 
autonomous  ground  vehicles.  Finally,  in  the  “Discussion”  chapter,  we  add 
intuition  to  our  results,  as  well  as  state  the  existence  of  two  interesting  bounds. 
Finally,  in  the  “Conclusion”  chapter,  we  conclude  with  recommendations  for 
future  research. 

Note,  in  the  baseline  model  developed  in  “Physical  Model  of  Transition,”  wave 
and  current  forces  are  initially  ignored.  As  no  model  currently  exists  for  the 
transition  that  can  facilitate  real-time  decision  making  by  an  autonomous  agent, 
the  main  focus  of  this  project  was  to  develop  a  simplified  model  for  initial  analysis 
that  in  the  future  could  be  expanded  and  refined  with  the  addition  of  a  more 
comprehensive  set  of  parameters. 
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Preliminaries 

Chapter  Overview 

In  this  chapter,  we  establish  terminology  and  develop  a  coordinate  system  that 
will  be  used  in  the  development  of  the  model  in  the  next  chapter.  We  follow  with 
a  precise  definition  of  the  problem  statement.  Then,  to  motivate  our  work,  we 
briefly  review  the  existing  literature  and  discuss  their  limitations  as  they  apply  to 
our  problem. 

Terminology 

We  define  the  map  m  c  R2  to  be  the  set  of  all  planar  points  where  the  vehicle's 
center  of  mass  could  be  positioned.  Define  a  terrain  elevation  function 
such  that  when  the  vehicle's  center  of  mass  is  at  w  e  mj  the  elevation  of  the 
ground  surface  or  sea  floor  at  w  is  E(w)  relative  to  the  sea-level. 

Define  a  draft  function  h  :  M  IR+  u  {0}  such  that  when  the  vehicle's  center  of 
mass  is  at  w,  the  draft,  or  vertical  distance  from  the  waterline  to  the  bottom  of  the 
hull,  is  /?(w).  Let  H  be  the  maximum  draft  of  the  vehicle  (corresponding  to  a 
floating  vehicle  at  sea). 

Now  define  the  world  '  to  be  the  surface  in  R3  that  denotes  the  position  of  the 
bottom  of  the  vehicle's  wheels/tracks.  More  precisely, 


W;=  {(J ■  3/’ -J  *  ^  (z,!0  fM.  :  =  max{£((a:,y)),-iy}  if  E((x,y))  <  0,  z  =  E((x,y))  if  E{(x,y))  >  0}. 


The  operating  environments  are  now  precisely  defined.  We  define  the  sea 
s  :=  e  w  I  =  if};  the  land,  «  -  ft***)  e  w  j  *«!,*))  =  o};  and  the 

transition  7* 1  *=  w  A((*,y))  e  ■ },  For  our  analysis,  we  consider  the 

closure  of  To,  denoted  ^  :=  Next  the  borders  are  defined. 

The  engage  line,  5nI-  where  S  is  the  closure  of  S.  Likewise,  the 
shoreline,  h  :=  ^  r  7  .  See  Figure  3  for  an  illustration  of  the  engage  line  and 
shoreline  relative  to  the  transition.  Then,  our  first  baseline  assumption  is  made: 

Assumption  1.  The  engage  line,  le,  and  shoreline,  /s,  are  straight  lines  and 
parallel  to  each  other. 

Without  loss  of  generality,  we  orient  the  coordinate  system  such  that  ~ 

g  V'  |  y  —  U}.  Qefjne  distance  between  the  engage  line  and  the 
shoreline: 
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y  :=  min  ||fo»lb)  ~  (^1,  »0ll 

^ II  ^ y  1 : -  I } G I c -i 

where  is  the  Euclidean  norm.  The  direction  from  the  engage  line  to  the 
shoreline  is  defined  to  be  the  positive  y-axis  direction.  That  is, 

h  =  {{r,y,2)  6  W  |  y  =  F} 

Another  baseline  assumption  is  made: 

Assumption  2.  The  sea-floor  is  planar  in  the  closure  of  the  transition.  That  is,  for 

.4  :=  {{z.y)  e  M  |  {i,  y,  €  T},  E  :  .4  -4  R  is  affine,  i.e., 

3  >  0  st  V  (x, y,z)  G  T and  V  yt  £  [0 ,F),  y2  £  (yi,F], 

Define  a  G  [(l  f  •  to  be  the  angle  of  inclination  of  the  sea-floor  in  the  transition. 

That  is,'1  :=  t;m  ,T  .  See  Figure  3  and  Figure  4  for  an  illustration  of  the 
transition  with  Assumptions  1  and  2. 


Figure  3:  Sideview  of  transition  with  Assumptions  1  and  2 


Figure  4:  Isometric  view  of  transition  with  Assumptions  1  and  2 
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Established  Coordinate  System 

In  this  section,  a  precise  coordinate  system  is  defined  for  the  analysis  that  will 
take  place  later.  Recall  that  a  point  in  T  is  an  ordered  triplet.  In  the  previous 
section,  we  implicitly  considered  a  coordinate  system  {x,y,z}  as  illustrated  in 
Figure  5.  For  example,  the  elevation  function  maps  a  point  in  the  {x,y}  plane  to 
the  z-axis.  This  coordinate  system  was  only  used  to  allow  for  intuitive  definitions 
of  the  vehicle's  operating  environments.  We  now  introduce  the  transformed 
coordinate  system  that  will  be  used  for  the  remainder  of  the  paper. 

Following  from  Assumption  2,  the  vehicle  maneuvers  in  the  plane  that  intersects 
the  shoreline  and  the  engage  line  as  depicted  in  Figure  3  and  Figure  4.  To  allow 
most  of  our  analysis  to  be  performed  in  two  dimensions  rather  than  three,  we 
define  a  new  coordinate  system  which  is  the  {x,y,z}  coordinate  system 

rotated  about  the  shoreline  by  2  '  ,  and  translated  in  the  direction  such  that 

X 

the  -axis  corresponds  to  the  engage  line,  as  shown  in  Figure  6.  Define  the  unit 
vectors  J,  and  k  corresponding  to  the  ^  V,  and  s  axis,  respectively. 


We  define  the  sea-floor  plane,  denoted  1  ,  as  the  closure  of  the  transition  rotated 
about  the  shoreline  as  is  done  for  the  {**$**}  coordinate  system.  That  is, 
f  :=  {(i.y.O)  €  B3  |  (i.ioG.o.0)  e  rt.  Similarly,  define 
distance  between  the  engage  line  and  shoreline  with  respect  to  the  sea-floor 
plane.  From  this  point  on,  we  always  refer  to  the  coordinate  system. 


Y  •=  ^ 

cosd  to  obtain  the 


Figure  5:  Position  coordinate  system 
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In  this  analysis,  we  consider  the  vehicle’s  heading  angle,  denoted  0,  which  is 

measured  relative  to  the  positive  *  -axis.  As  illustrated  in  Figure  7,  a  heading 
angle  of  0  corresponds  to  traveling  parallel  to  the  engage  line  and  the  shoreline, 
and  a  heading  angle  of  n/2  radians  corresponds  to  traveling  perpendicular 
towards  the  shoreline. 


Figure  7:  Heading  angle  6  with  respect  to  sea-floor  plane 


A  final  remark  about  notation.  A  symbol  denoted  with  a  bold  letter,  such  as  a 

position-dependent  force  F,  always  represents  a  vector  in  (\.e.,F :  T  r3), 
and  a  symbol  not  in  bold  which  has  a  corresponding  bold  symbol,  such  as  F, 


7 


Naval  Surface  Warfare  Center  Carderock  Division 
Center  for  Innovation  in  Ship  Design 
DUKW-21  Autonomous  Navigation:  Transitioning  Between  Sea  and  Land 

represents  the  magnitude  of  the  corresponding  bold  faced  force.  That  is  F(w)  := 
||F(w)||,  Vwe*\ 

Problem  Statement 

Before  presenting  the  problem  statement,  we  first  provide  some  context.  When  a 
wheel  or  tracked  amphibious  vehicle  makes  contact  with  the  engage  line,  a 
significant  amount  of  energy  is  lost  in  the  system.  Hence,  the  vehicle  effectively 
starts  from  near  rest  when  it  begins  its  journey  to  the  shoreline  from  the  engage 
line. 

Furthermore,  as  discussed  earlier,  the  tire/track  to  ground  surface  friction 
coefficient  is  very  small,  which  would  result  in  significant  tire/track  slipping, 
making  maneuvering  a  major  challenge.  It  may  be  most  practical  to  implement  a 
traction  control  system  on  the  vehicle,  which  would  keep  the  tires/tracks  rotating 
at  the  speed  that  would  maximize  the  tractive  force  generated  by  the  vehicle.  In 
this  case,  linear  paths  in  the  transition  are  preferred.  Now  a  formal  statement  of 
the  problem  is  given,  with  an  accompanying  illustration  in  Figure  8. 


Problem  statement:  Let  the  vehicle  start  at  position  1  '!1,  ^  ^  with  heading 

angle  w  *=  (^^7r),  initially  at  rest.  Let  ^  be  the  set  of  all  line  segments,  or  linear 
paths,  connecting  (0,0,0)  to  points  on /s.  That  is,  7  ^  €  L,  3  (A‘i ,  yi ,  0)  E 
Im,  Aeffi  sty  x  e  [CL  Til,  f;r,  Ai,0)  e  /.  Denote  the  time  it  takes  to  travel 

from  (0,0,0)  to  4  along  a  linear  path '  ^  ^ ■'  Our  objective  is  to  find 

argmin{t(/)} . 
idt 


Figure  8:  Set  of  linear  trajectories  from  p  on  the  engage  line  to  the  shoreline. 
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Related  Work 

To  our  knowledge,  there  has  been  no  research  on  path  planning  for  an 
amphibious  vehicle  in  the  transition.  In  particular,  there  is  no  existing  physical 
model  of  the  transition  that  has  low  enough  complexity  to  facilitate  the 
development  of  real-time  path-planning  algorithms  for  an  autonomous  agent  to 
use.  A  widely  used  physical  model  for  a  general  terrain  on  dry  land  is  developed 
in  (Rowe  &  Ross,  1990),  and  has  been  used  to  develop  successful  algorithms  for 
autonomous  navigation  on  dry  land  (Sun  &  Reif,  2005).  In  addition  to  not 
considering  buoyancy,  the  models  also  neglect  acceleration,  and  their  objective 
is  to  minimize  energy,  whereas  in  the  following  analysis  our  objective  is  to 
minimize  travel-time. 
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Physical  Model  of  Transition 

Forces  Generated  by  and  Acting  on  the  Vehicle 

In  this  section,  we  discuss  the  forces  in  the  system  consisting  of  the  vehicle 
driving  up  the  sea-floor  in  the  transition.  Recall  that  wave  and  current  forces  are 
neglected  in  this  baseline  model.  Note,  we  assume  that  all  forces  are  acting  on 
the  vehicle's  center  of  mass. 

The  only  considered  force  generated  by  the  vehicle  is  the  tractive  force  of  the 

ground  propulsors,  denoted  Ftrac  '■  T  — ¥  R3 .  We  assume  that  when  the 
vehicle’s  wheels  or  tracks  are  engaged  in  the  sea-floor,  the  efficiency  of  the 
waterborne  propulsors  will  be  negligible  compared  to  the  ground  propulsion  force 
and  is  omitted  in  our  analysis.  As  discussed  in  (Gillespie,  1992),  the  principal 
forces  acting  on  a  vehicle  driving  up  a  sloped  surface  on  land  are  the  weight,  the 
normal  force,  the  rolling  resistance,  and  the  drag,  denoted  W,  N.  Rrr,  D  :  T  K3, 
respectively.  In  the  transition,  there  is  an  additional  force  of  buoyancy,  denoted 
b  :  t  i3.  a  side-view  of  the  force  diagram  in  the  transition  is  shown  in  Figure  9. 


Figure  9:  Forces  generated  by,  or  acting  on  the  vehicle  in  the  transition 


Before  we  begin  the  analysis  of  the  transition  forces,  we  establish  consistent 
notation  to  precisely  describe  the  vehicle's  effective  weight.  We  define  the  weight 
W  :=  mg,  where  m  is  the  vehicle's  mass  and  g  is  the  gravitational  acceleration. 
Define  the  effective  weight,  We  :=W  +  B  (where  We  =  W  -  B).  The  effective  weight 
can  be  split  into  the  effective  normal  weight,  denoted  Wen,  and  the  effective 
tangential  weight,  denoted  Wet,  where  Wen  =  We  cosa  and  Wet  =  We  sina,  as 
shown  in  Figure  10.  Note  that  N  +  Wen=  0. 

It  then  follows  that  Rrr  =  crrWen,  where  crris  the  rolling  resistance  coefficient 
between  the  tires/tracks  and  the  ground  surface  (Gillespie,  1992). 
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Figure  10:  Effective  weight  (We)  split  into  normal  (Wen)  and  tangential  (Wet) 


Vehicle  Traveling  Perpendicular  to  the  Shoreline 

In  this  section,  we  consider  in  more  detail,  the  forces  generated  by,  and  acting 
on,  the  vehicle  as  it  drives  up  the  sea-floor  plane  in  a  straight  line  perpendicular 
to  the  shoreline  and  engage  line  (that  is,  0  =  tc/2).  A  more  general  discussion  for 
arbitrary  heading  is  presented  in  the  next  section.  The  following  assumption  is 
made: 

Assumption  3.  The  net  force  is  always  positive  in  the  positive  ■'  direction.  That 
is,  the  vehicle  can  always  overcome  resistance  and  opposing  forces  to  move 
forward/towards  the  shore. 

The  net  force  of  a  ground  vehicle  driving  perpendicularly  up  a  hill  with  an  incline 
of  angle  a  is  usually  modeled  with  Equations  (1)-  (3)  (Gillespie,  1992).  See 
Figure  1 1  for  an  illustration  of  the  sea-floor  planar  forces  in  the  transition. 


F  —  Ftrac  +  We<  +  Rrr  +  D, 


(1) 


F  =  Ftrac  ~  Wet  —  Rrr  —  D, 


(2) 


F  =  Fj 


(3) 
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Figure  11:  Sea-floor  planar  forces  in  the  transition 


Since  our  goal  is  to  find  the  linear  path  that  minimizes  travel  time,  without  loss  of 
optimality  we  assume  that  the  tractive  force  generated  by  the  vehicle  is  equal  to 
the  maximum  tractive  force  it  can  generate.  Hence,  for  the  remainder  of  the 
paper,  we  will  replace  FtracWith  Fmti-,  which  is  defined  as  the  maximum  tractive 
force. 

We  now  review  a  commonly  used  model  for  obtaining  the  maximum  tractive 
force  a  ground  vehicle  can  generate,  as  described  in  (Gillespie,  1992).  In 
general,  there  are  two  classes  of  properties  that  can  limit  the  maximum  tractive 
force.  The  first  is  engine  and  drivetrain  characteristics,  including  engine  torque, 
gear  ratios,  efficiency  of  the  drive  system,  and  the  inertia  of  the  engine  and 
drivetrain  components.  When  this  class  limits  the  maximum  tractive  force,  the 
expression  is: 


t'MTrl  =  - - - 


((/e  4-  4-  IdNf  4-  /u) 


where: 

Fmtm  is  the  maximum  tractive  force  the  vehicle  can  generate  when  the  upper 
bound  is  caused  by  the  first  class  of  properties 

Teis  the  maximum  engine  torque 

/Wis  the  combined  ratio  of  transmission  and  final  drive 

r/tf  is  the  combined  efficiency  of  transmission  and  final  drive 

r  is  the  radius  of  the  wheels 

/e  is  the  engine  rotational  inertia 

It  is  the  rotational  inertia  of  the  transmission  as  seen  from  the  engine  side 
ld  is  the  rotational  inertia  of  the  driveshaft 


(4) 
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Nf  \s  the  numerical  ratio  of  the  final  drive 

lw  is  the  rotational  inertia  of  the  wheels  and  axle  shafts 

a  is  the  acceleration  of  the  vehicle 

Equation  (4)  accurately  characterizes  the  maximum  tractive  force  for  large  friction 
coefficients.  For  small  friction  coefficients  between  the  tire  and  the  road,  the 
maximum  tractive  force  is: 


FmTtQ.  =  t'-fWe,  (5) 

where  pf  is  the  peak  friction  coefficient.  We  now  demonstrate  the  similarities  in 
the  structures  of  FMTri  and  FMTr2-  Introduce  the  constants: 


hi  ■  Tehtfritf 

r 

,  (Je  +  +  Id.Nf  +  Iw- 

h  :=  - * - 

Rewrite  Equations  (4)  and  (5)  as  follows: 


FMTri  —  77 y-  *  We  ~ 

vVe 

FmTt2  =  P-fWe  —  Oa  (6) 

Then,  the  general  maximum  tractive  force  the  vehicle  can  generate  can  be 
written  as: 

Fmtt  =  l*We  —  7  a,  (7) 

where: 

(I  =  f  ifc  FMTri  <  FmTt2 

\  i-tf  FMTri  >  FmTt2 

and: 

„  =  f  b-2  FMTri  <  F MTr2 
\  0  FMTri  >  F MTr2 

An  assumption  is  stated: 

Assumption  4.  The  track/sea-floor  peak  friction  coefficient  pf  and  the  track/sea¬ 
floor  rolling  resistance  coefficient  crrare  constant  throughout  the  transition.  It 
follows  that  p  is  constant  throughout  the  transition. 
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When  the  vehicle  travels  towards  the  shoreline  from  the  engage  line,  and  the 
heading  is  perpendicular  to  the  shoreline,  it  follows: 

F Af t v  -  (pWs  -  "ra)j  =  (p(W  -  B)  -  7a)  j  =  \pt(mg  -  pgVdup )  -  7a]j,  (8) 

where  FMTr  is  the  maximum  tractive  force,  p  is  the  fluid  density,  and  VdiSp  is  the 
volume  of  displaced  fluid. 

Moreover: 


Rrr  —  enj  —  ^tt  (ff  jBjcOSftj 

-  -  Crr  cos  a  (mg  -  pgVdiSp )  j . 


(9) 


Integrating  drag  into  the  model  significantly  increases  the  complexity.  Because  of 
the  expected  slow  speeds  that  the  vehicle  can  travel  at  in  the  transition,  it  was 
decided  that  drag,  which  is  proportional  to  the  speed  squared,  would  be  small 
relative  to  other  forces,  and  hence,  can  be  neglected  in  the  initial  analysis.  In  the 
future,  we  plan  on  studying  ways  of  integrating  drag  into  our  model  and  relaxing 
the  following  assumption: 

Assumption  5.  Drag  is  negligible  in  the  transition. 

We  couple  the  maximum  tractive  force  with  the  rolling  resistance  by  defining  the 
maximum  relative  ground  propulsive  force,  :=  Fjwrr  +Rrr  =*  Fmg  - 
Fmtt  -  Rrr-  Substituting  in  Equations  (8)  and  (9)  yields: 

FjtfG  =[(/*-  Crr  COS  Or)  (mg  -  pgVdiap)  -  7°]j  •  C 1[)) 

In  general,  precisely  calculating  the  volume  of  displaced  water  when  the  vehicle 
surfaces  is  very  challenging,  and  requires  intricate  knowledge  of  the  hull  form. 
Assuming  that  the  volume  of  displaced  fluid  is  a  linear  function  of  the  depth  of  the 
hull  greatly  decreases  the  complexity  of  the  problem,  so  for  now  we  make 
Assumption  6.  In  the  future  we  will  study  ways  of  approximating  the  hull  form 
piecewise  linearly. 

Assumption  6.  The  volume  of  displaced  fluid  is  a  linear  function  of  the  draft,  or 
vertical  distance  from  the  waterline  to  the  bottom  of  the  hull. 


Let  VHbe  the  maximum  displaced  volume,  corresponding  to  when  the  draft 

h.  Vt 

equals  H.  By  Assumption  6,  the  displaced  volume  at  a  draft  h  is  equal  to  ■,| 
We  can  now  extend  Equation  (10)  to  express  it  as  a  function  of  the  draft: 


H- 
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FMc{k)  =  (ft  - crr  cos  a)mg  -  7 a  -(ft-  crr 


cos  a)gp—h. 
±1 


(11) 


In  order  to  model  the  propulsive  forces  as  a  function  of  the  vehicles  position  in 

the  sea-floor  plane,  we  replace  the  draft  variable  h  with  the  variable 
corresponding  to  sea-floor  planar  distance  from  the  engage  line.  As  can  be  seen 

in  Figure  12,  ,,,,  =  ^  ~  y)  sina>  so  Equation  (11)  can  be  replaced  with: 


Fmg(v)  = 


Vjj  .  - 

(ft  —  crr  cos  a)mg  —  70  —  (ft  —  crr  cosa)gp-=-  sin(cv)(V'  —  y) 

H 


V, 


r 

=  (ft  —  crr  cos  a)  mg  —  (ft  —  crr  cos  a)gp^=r-  sin(cv)Y  —  7a 

XI 


Vj 


-1-0  -  cTT  cos  a) 9P^~  sin(a)y 


(12) 


Introducing  the  constants 

cti  :=  (p  —  crr  cos  o:)mg  —  (ft  —  crr  cos  a)gp^-  sin(a)V'  an(j 

a-i  \=  (p  -  Crr  cos  a) gp-£  sin  at  we  deliver  a  more  elegant  form  of  Equation  (12), 


Fmg(v)  =  ai  +  a2y  -  y a  (13) 

^  [1  0  FjV;f  g  ( y }  g  ( y)i  ■ 


To  review,  Equation  (13)  characterizes  the  sum  of  Ftrac(=  FMTr)and  Rrras  in 
Equation  (1).  Based  on  Assumption  5,  we  omit  the  drag  force  (D)  from  Equation 
(1).  Lastly,  we  integrate  the  effective  tangential  weight.  Recall, 
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tV  et(h)  =  We  sin  a  =  (TV  —  B)  sin«  =  (rug  —  pg^-h)  sin  a  (14) 

H 


Making  the  ^  ^  s^n  a  substitution  delivers, 


TVef(y)  =pi  +p2y ,  pi,p2  Gffi, 


(1G) 


where 


To  combine  all  the  discussions  together,  we  define  the  maximum  net  force  when 
the  vehicle's  heading  is  perpendicular  to  the  shoreline,  Fmax  :=  FMg  +  Wet,  where 
F max  —  FMg  -  Wet,  which  we  can  express  precisely  as  a  function  of  the  beach 


planar  distance  from  the  engage  line  (  ■'')  and  the  acceleration  (a): 


Fmax{if)  =  (oi  +  a2y)  -  (pi  +  p2y)  -  ja, 
where  =  Vma.T(y)j,  f0||0Wjng  from  Assumption  3. 


(17) 


Traveling  towards  Shoreline  with  Arbitrary  Heading 

Up  to  this  point,  we  have  only  considered  forces  as  the  vehicle  drives  towards 
the  shoreline  in  a  straight,  perpendicular  fashion.  We  now  add  to  our  analysis 
the  vehicle's  heading  angle,  0,  and  extend  our  discussion  to  the  vehicle  traveling 
at  a  non-normal  angle  relative  to  the  shoreline,  as  shown  in  Figure  12.  In  this 
section,  we  consider  the  forces  generated  by,  and  acting  on,  the  vehicle  when  it 
has  a  heading  angle  strictly  between  0  and  n.  In  other  words,  only  the  case  when 
the  vehicle  is  driving  towards  the  shoreline  is  considered  in  this  analysis. 

Equation  (17)  is  not  direction-dependent  because  it  follows  from  Equation  (1), 
which  assumes  the  vehicle  is  traveling  perpendicular  to  the  shoreline.  Another 
way  to  interpret  this  is  that  Wetwas  assumed  to  be  acting  on  the  vehicle  parallel 

jJ  JT 

to  the  direction  of  wheel/track  rotation.  When  ■  2  >  part  of  Wet  acts  on  the 
vehicle  perpendicular  to  the  direction  of  wheel/track  rotation,  which  has  a 
different  effect.  Wet  can  be  characterized  as  a  sum  of  a  vector  normal  to  the 
vehicle's  heading  direction,  denoted  Wetn,  and  a  vector  tangent  to  the  vehicle's 
heading  direction,  denoted  Wett,  both  also  contained  in  the  sea-surface  plane,  as 


16 


Naval  Surface  Warfare  Center  Carderock  Division 
Center  for  Innovation  in  Ship  Design 
DUKW-21  Autonomous  Navigation:  Transitioning  Between  Sea  and  Land 

illustrated  in  Figure  13.  Using  trigonometry,  it  can  be  shown  that  Wett=  Wetsin0 
and  Went  =  WetCOS0. 


Figure  13:  Effective  tangential  weight  split  into  vector  components  in  sea-floor  plane 

Went  causes  a  collinear  resistance  force,  Rsin  the  opposite  direction,  which  takes 
into  account  sideways  friction  of  the  tracks.  Tracks  typically  dig  into  the  soil, 
which  makes  it  even  harder  for  a  vehicle  to  slide  sideways.  Hence,  we  make  the 
following  assumption: 

Assumption  7.  The  vehicle  cannot  slide  sideways  in  the  transition  as  a  result 
of  gravity.  That  is,  Rs  =  -  Went-  See  Figure  14  for  an  illustration  of  this 
assumption.  _ 


Figure  14:  Sea-floor  planar  force  diagram  with  Assumption  7 
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To  review,  we  replaced  Wetwith  Wettin  Equation  (1).  Applying  a  similar  analysis 
as  in  the  previous  section,  we  can  extend  Equation  (13)  to  get  a  characterization 

of  the  maximum  net  force  of  the  vehicle  positioned  at  a  point  ■ ,r-  V 1  -  ^  and 

traveling  at  an  arbitrary  angle  ^  ^  Jl]- 


Fmaxix,  y,0)  =  (ai  +  a2y)  -  (pi  +P2y)sm0  -  7 a, 


(18) 


where  Fmax  (x ,  y,  0)  -  Fmax  (x,  y,  0)  (cos  i9i  +  sin  0j) . 


J7<  _ 

To  solve  for  the  maximum  acceleration,  set  ^maa?  —  ma 


a(x,y, 9}  = 


(ai  +  02$)  -  (pi  +  P2V)  sin  9 
m  +  7 


(19) 


di  '=  fH  do  =  ^  m  =  !>1  do  =  JJ2 

Defining  the  constants:  ’  !  '  ^  ™+r  ^  z  J 

we  get  a  characterization  of  the  maximum  acceleration  of  the  vehicle  positioned 
at  traveling  at  an  arbitrary  heading  angle  0: 


y,  9)  =  (di  +  a2y)  -  (pi  +  p2y)  sin  9 


(20) 
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Implications/Analysis  of  the  Force  Model 

Chapter  Overview 

In  this  chapter,  we  first  consider  the  problem  of  finding  the  linear  path  connecting 
an  initial  starting  position  on  the  engage  line  to  the  shoreline  which  minimizes 
travel  time.  This  is  an  important  problem  to  solve  because  the  amphibious 
vehicle  will  be  transporting  supplies  that  are  urgently  needed  ashore,  such  as 
food  and  medicine.  Because  of  the  low  friction  coefficient  between  the 
wheels/tracks  and  wet  sand,  the  maximum  propulsive  force  of  the  vehicle  is 
significantly  smaller  in  the  transition  than  on  dry  land.  Hence,  traveling  on  a 
slower  path  has  the  potential  to  greatly  increase  the  time  needed  to  deliver 
supplies. 

The  results  obtained  for  path  planning  in  the  transition,  as  we  show,  can  be 
applied  to  hills  on  dry  land  in  certain  situations,  which  we  discuss  later  in  the 
chapter. 

Note:  all  the  proofs  are  omitted  here  and  can  be  found  in  Appendix  A. 


Fastest  Linear  Path  in  Transition 

In  this  section,  we  present  the  main  contribution  of  this  report,  which  is  a  closed 
form  expression  for  the  heading  angle  corresponding  to  the  fastest  linear  path  for 
a  vehicle  to  travel  from  an  initial  starting  point  in  the  engage  line  that  terminates 
at  the  shoreline.  The  expressions  are  stated  in  Theorems  5  and  9.  From 
Assumption  7,  we  get  Equation  (20),  which  demonstrates  how  the  maximum 
acceleration  the  vehicle  can  obtain  is  dependent  on  the  heading  angle  of  the 
traveling  vehicle. 

The  critical  step  in  finding  the  optimal  heading  angle  is  Lemma  1 ,  which  yields  a 
closed  form  expression  for  the  travel  time  along  any  linear  path  assuming  the 
vehicle  travels  as  fast  as  possible.  The  key  to  proving  this  lemma  is  to  recognize 

A  A 

ai  _  pi 

that  which  we  establish  below. 


Lemma  1.  Let  the  vehicle  start  at  position  (0,0,0)  e  4  with  heading  angle 
6  e  (0,7r) ,  initially  at  rest.  Let  /  be  a  straight  line  segment  along  0  connecting 


(0,0,0)  to  a  point 


(Ai.r.o) 


€  l 


where  taLi* 


>0. 


Suppose  Assumptions  1 
through  7  hold,  and  that  the  vehicle  will  always  travel  as  fast  as  possible.  Also, 

assume  that  neither  n 2  nor  equals  0.  Then,  the  travel  time  along  line  /  is: 
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t  (  A)  =  K 


1  +  i? 


Where  A  =  =  tan/9,  K  = 


^  ft  2  yj  1  +  —  f>2 


dy 


„  0,  and  k  =  =  U- 

V(S+*)2-i;s 


(21) 


The  last  hypothesis  in  Lemma  1,  that  neither  ■  '■  nor  -  equals  0,  implies  that 
the  displaced  water  volume  is  nonzero,  which  means  the  vehicle  is  in  the 

transition  and  not  on  dry  land.  Observe  that  if  "  “  =  0,  Equation  (21)  does  not 
make  sense.  We  could  have  weakened  the  hypothesis  to  just  say 
however  «2  =  0  =*  =  0  . 


£2  >  1 

Theorem  2.  Let  the  conditions  in  Lemma  1  hold.  Further,  assume  that  «2  2  • 

Then,  the  linear  path  which  minimizes  travel  time  corresponds  to  the  vehicle 
traveling  at  angle: 


=  tan  1 


\  4  ( — «a<s — Y  _  i 

\  \  \  H—Crr  COSO  ) 


(22) 


To  clarify  the  last  hypothesis,  observe  that  if  ^2  ft  —  cTT  cos  oe  -  b  Equation  (32) 

does  not  make  sense.  Intuitively,  this  means  that  Theorem  2  is  applicable  when 
the  friction  coefficient  between  the  tires/tracks  and  ground  surface  is  sufficiently 
small  and  the  rolling  resistance  coefficient  is  sufficiently  large  (such  as  on  wet 
sand),  and/or  the  slope  of  the  sea-floor  plane  is  sufficiently  large. 


We  now  present  an  example  of  a  situation  when  Theorem  2  applies,  that  is, 

h.  >  i 

when  “a  2’  Heavy  trucks  on  sand  generally  encounter  a  rolling  resistance 
coefficient  of  about  0.25,  so  assume  crr  =  0.25  (Gillespie,  1992).  Further,  wet  tires 
can  have  friction  coefficients  of  approximately  0.4,  so  assume  p  =  0.4  (Jones  & 
Childers,  2000).  Now  observe  that  when  the  sea-floor  plane  has  an  angle  of  a  = 

P2  _  sin  Q  _  sill  10  ^  -|  -|  O  ^  l_ 

^  qo  a 2  crr  cosa  0.4— 0.25  cos  10  ^  *  2* 

Then,  by  Theorem  2,  the  optimal  heading  angle  p*  ~  >  which  is  noticeably 

different  than  the  perpendicular  heading  angle  of  90°. 

h.  <  k 

Now  we  consider  the  case  not  covered  by  Theorem  2,  that  is,  when  —  2  ■ 
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Theorem  3.  Let  the  assumptions  in  Lemma  1  hold.  Further,  suppose  that 

A  -i 

£1  <  I 

_  2'  Then,  the  linear  path  which  minimizes  travel  time  is  the  normal  line 
connecting  (0,0,0)  to  ls.  That  is,  0t  =  f . 


Fastest  Linear  Path  on  Dry  Hill 

Observe  that  in  Lemma  1,  Theorem  2,  and  Theorem  3,  it  is  assumed  that 

A 

and  a'2  are  not  equal  to  zero.  This  makes  sense  when  the  vehicle  is  in  the 
transition,  and  the  displaced  water  volume  is  strictly  greater  than  zero.  However, 
as  we  show  in  this  section,  the  final  characterization  of  the  optimal  linear  path  in 
the  transition  does  not  depend  on  the  effective  weight.  In  other  words,  the 
buoyancy  term  cancels  out,  and  our  results  can  be  extended  to  traveling  on  dry 
land.  We  demonstrate  this  in  the  remainder  of  the  chapter. 

Consider  a  vehicle  travelling  up  a  sloped  terrain.  Similarly  to  the  definition  of  the_ 
closure  of  the  transition,  T,  we  define  a  planar  surface  of  a  hill,  denoted  byH.H 
is  assumed  to  be  an  affine  two-dimensional  surface  with  a  inclination  angle  a, 
and  contains  the  closure  border  lines:  the  bottom  line,  ^  c  and  the  top  line, 
h  c  -v  analogous  to  le  and  ls  for  T.  Assumptions  1  through  4  and  7  are  also 
extended  to  hold  for  H- 

Observe  that  for  a  planar  surface  of  a  hill,  our  discussion  and  results  from  the 
previous  two  chapter  still  hold,  with  the  exception  that  now  Vh  =  0,  which  implies 

that  “2  =  *2  =  P2  =  m  =  o, ,  and  that  a, =  (,i "  Crr  cos  0  )ms  and  pi  =  m»sin  ”■  The 

results  are  now  stated. 


Lemma  4.  Let  the  vehicle  start  at  position  (0,0,0)  e  lb ,  with  heading  angle 
9  e  (0,;r) ,  initially  at  rest.  Let  /  be  a  straight  line  segment  along  0  connecting 

(0,0,0)  to  a  point  €  lt'  where  -  ^  Suppose  that  the  vehicle 

will  always  travel  as  fast  as  possible.  Then,  the  travel  time  along  line  /  is: 


t{\)  s  sfV'F) 


1  + 


\  2*!  ^1  + 


JL 

~~n 

u 


(23) 


Theorem  5.  Let  the  assumptions  in  Lemma  4  hold.  Further  suppose  that 

£1  >  1 

-'I  2'  Then,  the  linear  path  which  minimizes  travel  time  corresponds  to  the 

vehicle  traveling  at  angle  <9*  as  characterized  in  Theorem  2,  Equation  (22). 
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Theorem  6.  Let  the  assumptions  in  Lemma  4  hold.  Further,  suppose  that 

A  .1 

EL  <  I 

a  i  —  2 '  Then,  the  linear  path  which  minimizes  travel  time  is  the  normal  line 
connecting  (0,0,0)  to  lt.  That  is,  <9*  =  f . 
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Discussion 


In  Theorems  2  and  3,  we  characterized  the  heading  angle  corresponding  to  the 
optimal  linear  path  from  a  point  in  the  engage  line  to  the  shoreline  when  the 
vehicle  starts  at  rest.  Additionally,  we  gave  a  realistic  scenario  when  the  optimal 
heading  angle  is  significantly  different  than  the  angle  perpendicular  to  the 
shoreline. 


The  characterization  of  the  optimal  heading  angle  in  Equation  (22)  has  an 
intuitive  interpretation.  Define  the  constant: 


Now  rewrite  Equation  (22)  as: 


sin  a 

p  —  crr  cos  a 


=  tan  1 


(23) 


(24) 


Observe  that  for  constant  p  and  crr,  as  a  which  makes  Equation  (24) 

>  ,  1  jj  y  T 

not  suitable.  However,  according  to  Theorem  2,  as  s  2’  -  •  Intuitively, 

this  can  be  understood  as  saying  that  for  fixed  p  and  crr,  there  exists  a  sufficiently 
small  a,  not  necessary  equal  to  0,  such  that  the  fastest  linear  path  is  the 
perpendicular  path. 

On  the  other  hand,  observe  that  for  constant  p  and  crr,  as  Q'  f 1  £  h  °- 
Intuitively,  this  can  be  understood  as  saying  that  for  fixed  p  and  crr,  there  exists  a 
strictly  positive  lower  bound  to  the  optimal  heading  angle  for  linear  paths, 
regardless  of  a.  In  other  words,  at  a  certain  point,  reducing  the  effect  of  effective 
tangential  weight  will  not  save  time. 

In  “Implications/Analysis  of  the  Force  Model”,  it  was  shown  that  the  optimal  path 
had  the  same  form  whether  in  the  transition  or  on  a  hill  on  dry  land.  In  other 
words,  the  effective  weight  has  no  effect  on  the  orientation  of  the  optimal  linear 
path.  The  key  to  understanding  when  the  optimal  linear  path  will  be  perpendicular 

c  i 

to  the  shoreline  or  not  is  the  C,  term.  When  s  2>  the  fastest  linear  path  is  not 
perpendicular.  This  happens  when  the  incline  angle  is  sufficiently  large,  and 
more  importantly,  the  difference  between  the  friction  coefficient  and  the  rolling 
resistance  coefficient  is  sufficiently  small.  This  can  happen  in  the  transition,  as 
well  as  on  sloped  land  that  has  a  slippery  surface. 
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Conclusions 

Summary 

Several  accomplishments  were  achieved  in  this  report.  An  original  baseline 
model  was  proposed  for  the  transition,  which  is  the  operating  region  where  an 
amphibious  vehicle's  wheels  or  tracks  are  engaged  with  the  sea-floor,  at  the 
same  time  as  the  vehicle  itself  is  partially  submerged  in  water.  The  model 
developed  takes  into  account  the  ground  propulsive  force,  gravity,  buoyancy, 
sliding  friction  between  tires/tracks  and  sea-floor,  rolling  resistance  of  tires/tracks 
and  sea-floor,  and  engine  and  drivetrain  characteristics,  including  engine  torque, 
gear  ratios,  efficiency  of  the  drive  system,  and  the  inertia  of  the  engine  and 
drivetrain  components.  The  model  ignores  drag,  waves,  currents,  and  other  surf 
zone  forces,  which  are  to  be  addressed  in  future  work.  It  also  assumes  ideal 
properties  about  the  shape  of  the  transition. 

Using  the  model,  novel  observations  were  made  about  route  optimality.  Consider 
C,  =  sina/(p  -  crrcosa),  where  a  is  the  angle  of  inclination  of  the  transition  hill,  p  is 
the  peak  sliding  friction  coefficient  between  the  tires/tracks  and  the  sea-floor,  and 
crr  is  the  rolling  resistance  coefficient  between  the  tires/tracks  and  sea-floor.  It 
was  shown  that  when  C,  >  14,  the  fastest  path  for  the  vehicle  to  traverse  in  the 
transition  is  not  perpendicular  to  the  shoreline.  Moreover,  realistic  cases  were 
found  when  C,  >  14.  For  this  case,  a  closed-form  characterization  for  the  heading 

angle  corresponding  to  fastest  linear  path  was  found:  =  ,  where  a 

heading  angle  of  0°  corresponds  to  traveling  parallel  to  the  shoreline  and  a 
heading  angle  of  90°  corresponds  to  traveling  perpendicular  towards  the 
shoreline. 


Future  Work 

In  our  current  model,  we  do  not  impose  an  upper  bound  for  vehicle  speed.  In 
reality,  in  addition  to  the  position  and  heading  angle,  the  maximum  acceleration 
is  also  dependent  on  the  vehicle's  speed,  so  Equation  (20)  does  not  fully  capture 
the  maximum  acceleration.  This  factor  should  be  considered  in  the  future. 

In  calculating  the  travel  time  along  a  linear  path,  it  was  assumed  that  the 
vehicle's  initial  speed  was  zero.  In  the  future,  we  plan  to  relax  this  assumption 
and  evaluate  the  optimal  linear  path  for  a  general  starting  speed.  Then,  even  if 
the  starting  speed  is  zero  at  the  engage  line,  as  the  vehicle  travels  up  the 
transition,  the  optimal  linear  path  from  its  present  position  to  the  shoreline  will 
change.  Even  though,  as  stated  earlier,  maneuvering  precisely  in  the  transition 
would  be  very  challenging,  a  piecewise  linear  path  could  still  be  feasible. 
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Besides  loss  of  optimality,  there  is  a  potential  loss  of  feasibility  with  a  strictly 
linear  path.  Throughout  this  paper,  we  assumed  that  there  is  no  bound  on  the 
;E-axis  of  feasible  "land  entry  points."  If  the  bound  were  present,  Theorems  5 
and  9  would  yield  unfeasible  trajectories.  We  plan  to  address  this  limitation  in 
future  research. 

More  work  must  be  done  to  add  realism  to  the  model  proposed  in  this  report.  In 
particular,  wave  and  current  forces  in  the  transition  are  likely  to  noticeably 
influence  the  dynamics  of  the  vehicle.  Also,  the  transition  might  not  be  a  planar 
hill,  and  more  general  transition  geometries  should  be  incorporated  into  the 
model.  Similarly,  the  displaced  water  volume  is  not  necessarily  a  linear  function 
of  the  draft,  so  analysis  incorporating  specific  hull  forms  into  the  model  to 
calculate  the  buoyancy  more  precisely  is  needed.  Finally,  it  is  ultimately  desired 
to  integrate  obstacles  into  the  model  because  on  a  real  beach,  not  all  paths  will 
be  feasible. 


Recommendations 

The  model  and  algorithms  discussed  in  this  report  are  ready  to  be  tested  in 
controlled  environments.  The  model  beach  would  have  to  be  close  to  planar  and 
the  water  would  have  to  be  calm.  Before  testing  this  model  on  DUKW-ling, 
various  tests  need  be  done  on  the  DUKW-ling  itself  to  characterize  certain 
properties,  such  as  the  peak  friction  coefficient  and  rolling  resistance  between 
the  wheels  and  sea-floor,  which  are  pieces  of  information  that  are  taken  for 
granted  in  this  report. 
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Appendices 

Appendix  A  -  Proofs 

Proof  of  Lemma  1.  Before  we  begin  the  proof,  first  recall  that  for  a  position 
dependent  acceleration  along  a  straight  line  /,  denoted  afsj,  where  s  el,  the 
velocity  at  s  is  equal  to  (Pestel,  1968): 


■u(s)  = 


[  25 } 


where  v(s)  is  the  velocity  at  position  s,  and  v(so)  is  the  velocity  at  some  initial 
position  so.  The  approach  we  use  in  this  proof  is  to  interpret  the  integral  term  in 
Equation  (A1)  as  a  line  integral,  and  then  to  find  convenient  bijective  (one-to-one 
and  onto)  parameterizations  to  simplify  the  calculations. 


Indeed,  since  /  is  a  straight  line  segment  in  the  plane  spanned  by  the  ^  and 
'  axis,  we  know  that  the  following  is  a  bijective  parameterization  of  /: 
r  :  [0.  A.i]  ->  I defined  by  r  (£}  =  (5,  Ai),  where  \  ” 

seen  in  Figure  15, 


As  can  be 


sij  i(f  — 


I 

V 


+  Y2 


v 


n  + 


(26) 
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Figure  15:  Arbitrary  linear  trajectory,  (/),  from  origin  to  point  on  shoreline. 


Substituting  Equation  (26)  into  Equation  (20)  yields  a  characterization  of  the 
acceleration  on  an  arbitrary  linear  path  with  a  corresponding 


-  Qi  +  a-2V  - 


Pi  +  M 


For  an  arbitrary  (u,v,E((u,v)))  el ,  let  l(uv)  c=  /  be  the  straight  line  segment 


(27) 


connecting  (0,0,0)  to  ( u,v,E((u,v ))).  Since,  by  assumption,  the  tractive  force  of  the 
vehicle  is  always  equal  to  its  maximum  tractive  force,  and  the  initial  speed  is  0, 
from  Equation  (25),  we  have: 


(it.  v.  \)  =  f  a{r)dT 

V 


(28) 


Observe  that  restricting  r  to  the  domain  u ]  ^  1  is  a  bijective 

.  .  ..  f(  M  ,  ..  i^(x)  =  ( I .  A)  =>  ||4^  (:r)  1 1  =Vl  +Aa. 

parameterization  of /(U,v).  Note  that  dx  ^  '  MdxWM 

By  definition  of  line  integral,  it  follows: 


r 

ry/* 

dr 

a(r)dr  =  f  a  (r  (£) ) 

Q'l 

dx, 


(29) 


which,  from  Equation  (28),  implies: 
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/  - - -  rSM 

v  (;r,  y,  A)  =  i  /  2\/ 1  4-  A2  /  a  (ir,  A;c,  A)  dx 


(30) 


Evaluating  the  integral  using  Equation  (27)  yields: 


Jo  Jo 


*/ 1  +  i 


-l 

-l 


0.1  - 


Pi 


-I-  I  do  A  — 


PsA 


\/1  +  ^ 

Pi 

\/  1  +  A^ 


+  do  A 


V  J  +  ^ . 

Pa  A 
V/ 1  +  AT 


21 


37 


fll  - 


fil  - 


Pi 


\ 


x  +  (  a2A  — 


■p2  A 


\/ 1  +  i 


:r 


<ii': 


dx 


I  y/A 

lo 


Pi 


\/  J  +  A*  ^  A 


I  Y  +  - 7====  I  To  (3D 


^  1  +  A* 


r  a2 


So: 


v(x,y,A)  = 


2VI  +  A2 
A 


fli 


Pi 


\A+a 


2?  +  «2 


P2 


\ 


I 


(32) 


4i-  =  |i,  *.•  =  Cg- - y 

The  critical  observation  is  that  a<1  p2  which  we  call  p3  n  B111  Q 
This  implies: 
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v  (x,y.  A) 


V  ky  +  y2 


( 33 ) 


Now  we  make  another  bijective  parameterization  of  /,  called 

defined  by ~  (M  -  * ®  -  (i- 0  -  W ® 

using  line  integrals,  we  now  solve  for  travel  time: 


io, y] 


L 


-  )f]  +  i 


Again, 
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Proof  of  Theorem  2.  We  differentiate  the  result  of  Lemma  1 ,  Equation  (21 ): 


dt 

dX 


(A)  = 


\ 


3* 


/- 


1+- 


(42/TT J-h)  (^)  -(!+£)  (-0-)  (j?) 


(a2  i^/l  +  x?  -  P2) 


-  (H 

-  (H 

-  (H 


(«2y"l  +  jfc-fe)  ' 

\/  1  +  5" 


1/2  2  |  P3-«2V  1  +  ^  \  +  iiafl  +  tW) 


A3 


A3v' 


(oa^/l  +  jpr-Pa) 


2p2  - 


202^1+^+02^1  +yr 

A3  \/  1  +  ^  («-2^/l  +  -  P2) 

2p2  -02^/1  +  jr 

(a2\/  1  +  jpr  -  P2) 


To  find  the  critical  value,  denoted  A» ,  we  set  j|(A*)  =  0 ,  which  implies: 


~P2  ~  &2  \  /  I  +  =  1 } 


.  +  -  =  2^ 
/  A^  «2 


1  +  A?  4 


(I) 


\ 


/ 


(35) 


a! 


A,  = 


(36) 


(37) 


Suppose  A*  were  a  minimizer  of  t.  then,  there  must  be  a  corresponding  ^  * 

such  that  Ax*  1 ' ' 1 '  'J  :H  ‘  where  <9,  is  the  heading  angle  corresponding 

to  the  linear  path  associated  with  A* ,  or  the  optimal  linear  path.  Equation  (22) 
then  immediately  follows,  which  is  what  we  wished  to  show. 


It  remains  to  be  shown  that  A*  is  indeed  a  minimizer  of  t.  It  suffices  to  show  that 
A  >  A*  =>  J/(A)  >  0  and  A  <  A*  =>  Jj-(A)  <  0 .  Since  the  denominator  in  Equation 
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(35)  is  strictly  positive,  it  suffices  to  show  that  A  >  A,  =>  '^2  V 1  +  x*  -  02 11  +  xO  >  0 
and  A  <  A,  =>  “/J2  V  1  +  **  02  +  ^  U‘  Indeed,  suppose  A  <  A,.  It  follows: 


Similarly,  for  A  >  A, : 


2p2 


(38) 


(39) 

□ 
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Proof  of  Theorem  3.  Observe: 


h 

«2 


2 


1  Pi  <  «2 


2 pi 


/ 


+  yry  ^  &2 


SP2 


<  II. 


(40) 


As  can  be  seen  in  Equation  (35),  the  denominator  of  jj(A)  is  strictly  positive,  so 
Equation  (40)  implies  that  Jj-(/l)  is  nonincreasing,  which  means  that  as  A  -» oo , 

f(A)  gets  smaller  or  stays  the  same,  which  means  the  line  with  slope  oo ,  the 
perpendicular  line,  minimizes  travel  time. 


□ 


Proof  of  Lemma  4.  Everything  in  the  proof  for  Lemma  1  holds  up  through 
Equation  (32).  Now  substitute  =  ^2  =  ^  to  get: 


v  (x,y,\)  = 


2^/TTA2 

A 


«i  - 


Pi 


1  + 


v[xty,X)  =  y/yd2ai^j 

Applying  a  process  similar  to  Equation  (34)  yields: 


+  j?-'2pi 


which  implies  Equation  (23). 


(41 ) 


(42) 


□ 
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Proof  of  Theorem  5.  Differentiating  the  result  of  Lemma  4,  Equation  (23),  yields: 


dt  _  V?  (2pi  —  &lyj\  +Jt) 

"  =  A3v^T(,>alv^T-#1)3/2 

To  find  the  critical  value,  denoted  X, ,  we  set  -^(X,)  =  0 ,  which  implies: 


) 


(43) 


2pi 


-  ai 


V 


+  a?  =l>=i' 


1 

a! 


A*  = 


(44) 


(45) 


Suppose  2*  were  a  minimizerof  t.  then,  there  must  be  a  corresponding  ^  '  * 

such  that  Ai*  l',n  r  where  <9, is  the  heading  angle  corresponding 

to  the  linear  path  associated  with  X,  or  the  optimal  linear  path.  Equation  (22) 
then  immediately  follows,  which  is  what  we  wished  to  show. 


It  remains  to  be  shown  that  X,  is  indeed  a  minimizer  of  t.  It  suffices  to  show  that 
X  >  X,  =>  jj(X)  >  0  and  X  <  X,  =>  jj(X)  <  0 .  Since  the  denominator  in  Equation 


(43)  is  strictly  positive,  it  suffices  to  show  that  X  >  X, 


-p‘2 


\J  1  +  J?  -  “2  (l  +  -&)  >  0 


and  X  <  X* 


V 


m 


(1  +  i)  <  o. 


Indeed,  suppose  X<X.  It  follows: 
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Similarly,  for  X  >  A* : 


\2  \2 
X-  >  A: 


I 


4(fr)  >^+1^2{fl)>\ll  +  ^ 

r 


2(tL)  ./1  +  ±>±  + 
‘  Va;,y  V  A-  "  A2 


V  1  +  A2  _  Al  l  1  +  A2  j  >  ° 


(47) 

P 


Proof  of  Theorem  6.  Using  a  similar  procedure  as  the  proof  in  Theorem  3,  one 
can  establish  Theorem  6. 

□ 
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